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Solubility of CaHP0 4 * 2H 2 in the Quaternary System 
Ca(OH) 2 - H 3 P0 4 - NaCI - H 2 at 25 °C * 
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Solubility of CaHP0 4 ■ 2H 2 was determined in the quaternary system Ca(0H) 2 — H 3 P0 4 — NaCI — 
H 2 at 25 °C in the pH range 4.39-6.38; ionic strengths of the saturated solutions varied from 0.00485 to 
0.545, Satisfactory constancy in the solubility product, (K sp = [Ca 2+ ] • [HPO* ] • y Ca a+ • yniw- = 
2.49 ± 0.05 X 10" 7 mol 2 Xl -2 ) was obtained when (i) the ion activity coefficients, y„ were calculated 
with the Debye-Huckel equation, log y,- =—y4Z?V7/(l + Ba,-V7) +0.0626 /; the value 0.0626 for th 
coefficient in the linear term was derived from the solubility data by utilizing a statistical procedure, a 
(ii) formation of an ion pair NaHPOj was taken into account; a statistically derived value for t 
stability constant of this ion pair is 7.0 ± 2.4 1 X mol 1 . The ion pair NaHP() 4 appears to have significant 
concentrations in physiological fluids. 
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1. Introduction 

The sparingly soluble calcium phosphates, 
Ca 5 (P0 4 ) 3 OH, Ca 8 H 2 (P0 4 )«-5H 2 0, /3-Ca 3 (P0 4 ) 2 , 
CaHP0 4 -2H 2 and CaHP0 4 are of considerable 
physiological interest. For this reason, we have been 
investigating their solubilities with two major objec- 
tives in mind. First, we wish to develop an internally 
consistent set of thermodynamic quantities — solu- 
bility products, enthalpies and entropies of dissolu- 
tion—for as many of these salts as is feasible experi- 
mentally. Second, we wish to try to develop ionic 
models which will satisfactorily explain the experi- 
mental solubility data, and which can subsequently 
be used to describe the behavior of these salts in 
physiological systems. 

These data were collected at four temperatures 
(5, 15, 25, 37 °C), over a broad range of concentration 
(approximately 1000-fold variations in H+, Ca, and P 
concentrations), and for four salts (Ca 5 (P0 4 ) 3 OH, 
£-Ca 3 (P0 4 ) 2 , CaHP0 4 -2H 2 and CaHP0 4 ). Nearly 
all of these measurements were carried out in the 
ternary system, Ca(OH) 2 -H 3 P0 4 — H 2 0, with the 
objective of using the simplest possible system to 
assure ease of interpretation of the data [1-6]. * In 
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only one instance was a study made in a quaternary 
system (Ca 5 (P0 4 ) 3 OH in the system Ca(OH) 2 -H 3 P0 4 - 
H 2 — HC1) [7]. In the present study, solubility 
measurements were extended to include the physio- 
logically important salt, NaCI, as the fourth component. 

In order to calculate the ion activity products a 
knowledge of ion activity coefficients is required; 
these are usually calculated with the use of either the 
Debye-Hiickel equation [8] , or the Davies equation 
[9] when the solutions have low ionic strengths. 
Previous solubility studies on calcium phosphates 
[1-7] have shown that both of these equations for 
calculating ion activity coefficients are adequate; the 
ionic strengths of most of these solutions were less 
than 0.1. 

The models which satisfactorily account for the 
above data require the assumption of the presence of 
the ion pairs CaHPOJ and CaH 2 PO+ [1, 2] in addition 
to the usual dissociation products of Ca(OH) 2 , H 3 P0 4 , 
and H 2 0. In the present work we find it necessary to 
include the ion pair NaHP0 4 . 

The purposes of the present investigation were 
(i) to test the usefulness of the expressions for ion 
activity coefficients in the system Ca(OH) 2 -H 3 P0 4 - 
NaCl-H 2 at ionic strengths considerably higher than 
0.1, and (ii) to determine whether or not there is an 
interaction between the ions of NaCI with those of 
Ca(OH> 2 and H 3 P0 4 . This was achieved by determina- 
tion of the solubility of CaHP0 4 • 2H 2 in the quater- 
nary system over a range of NaCI concentrations up to 
0.5 M. The salt CaHP0 4 • 2H 2 was selected for this 
purpose because (i) it is relatively easy to prepare in 
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a pure state, (ii) its relatively high solubility under the 
most propitious experimental conditions contributes 
to accuracy in the derived constants, and (iii) our 
experience indicates that its chemical properties are 
conducive to rapid attainment of equilibrium between 
it and the aqueous phase. 

2. Materials and Methods 
2.1. Materials 

The CaHP0 4 * 2H 2 used in this study was prepared 
by bubbling ammonia into an aqueous solution satura- 
ted with reagent grade Ca(H 2 P0 4 )2 * H 2 maintained 
at about 5 °C [10]; crystals that precipitated in the 
pH range 4.0 to 4.2 were collected, washed with 
distilled water and acetone, air dried, and stored in a 
refrigerator to prevent dehydration until needed. 

The H3PO4 solutions of various concentrations 
(Po, table 1) were made from doubly crystallized 
H3PO4 * iH 2 and freshly boiled distilled water. 
Known weights of reagent grade NaCl, dried for 
24 hours at 125 °C, were dissolved in the H3PO4 
solutions to provide the desired concentrations of 
NaCl. 

2.2. Solubility Experiments 

Each 2.5 g sample of CaHP0 4 ■ 2H 2 was equili- 
brated with 200 ml of NaCl-HjPC^ solution for two 
week periods in a thermostated water bath at 25 °C 
±0.01 °C equipped with an end-over-end rotating 
device. The equilibrated solutions were filtered 
through cellulose filters with average pore diameters 
of 0.24 fjLm. The pH values of the solutions were 
measured with a research pH meter which was stand- 
ardized with NBS certified buffers; the pH cell was 
{glass electrode | test solution || saturated KC1 1 Hg 2 Cl 2 ; 
Hg}; the estimated standard error in the pH measure- 
ments is ±0.01 pH units. Calcium concentrations, 
Ca*, of the equilibrated solutions were determined 
with a relative error of 0.025 by means of an atomic 
absorption spectrophotometer. Phosphorus concen- 
trations, P t , were determined to a relative error 0.015 
with a double beam spectrophotometer using vanado- 
molybdate reagent [11]. 

2.3. Calculations 

a. Ion-Activity Coefficients 

In most of our calculations, the Debye-Huckel ex- 
pression, eq (1), was used to evaluate the activity 
coefficients, yu needed to calculate ion activity 
products. 

\o gyi ={-AZ\ Vil(l+B-ai v7)}+0/. (1) 

In this equation A and B are temperature dependent 
constants; Zi is the charge number of the ith ion; a\ 
is its ion-size parameter and / is the ionic strength, 



Table 1. Solubility o/CaHP0 4 -2H 2 in the system 
Ca(OH) 2 -H 3 P04-NaCl-H 2 at 25 °C 


Solution 

No. 


NaCl, 
M 


Equilibrium Composition 
of Solutions 


Ca* 


pH 


Ca* 
MxlO 3 


P. 

MxlO 3 


Pt-Po 



: 14.8 MX 103 



A-l 


0.000 


4.39 


14.7 


29.8 


0.980 


A-2 


.050 


4.48 


14.8 


29.8 


.987 


A-3 


.100 


4„54 


14.9 


29.8 


.993 


A~4 


.150 


4.58 


15.0 


29.8 


1.00 


A-5 


.250 


4.64 


15.1 


29.8 


1.01 


A-6 


.500 


4.70 


15.2 


29.9 


1.01 



P =13.2MX10* 



B-l 


0.000 


4.47 


13.2 


26.2 


1.02 


B-4 


.150 


4.68 


13.5 


26.3 


1.03 


B-6 


.500 


4.80 


14.0 


27.2 


1.00 



P =10.5MX103 



C-l 


0.000 


4.62 


10.6 


21.8 


0.938 


C-2 


.050 


4.74 


10.7 


21.9 


.939 


C-3 


.100 


4.81 


10.8 


21.9 


.947 


C 4 


.150 


4.85 


10.9 . 


21.8 


.965 


C-5 


.250 


4.91 


11.1 


22.0 


.965 


C-6 


.500 


4.97 


11.4 


22.1 


.983 



Po-3.52 MX10 3 



D-l 


0.000 


5.41 


3.74 


7.30 


0.989 


D~2 


.050 


5.59 


3.93 


7.47 


.995 


D-3 


.100 


5.66 


4.08 


7.59 


1.00 


D-4 


.150 


5.70 


4.24 


7.76 


1.00 


D-5 


.250 


5.74 


4.52 


8.02 


LOO 


D-6 


.500 


5.79 


4.98 


8.74 


0.954 



Po = 1.21 MX10 3 



E~l 


0.000 


6.19 


1.60 


2.79 


1.01 


E-2 


050 


6,32 


2.06 


3.23 


1.02 


E--3 


.100 


6.36 


2.36 


3.51 


1.03 


E-4 


.150 


6.37 


2.57 


3.74 


1.02 


E~5 


.250 


6.38 


2.93 


4.12 


1.01 


E-6 


.500 


6.37 


3.62 


4.82 


1.00 



I=^CtZ% 



(2) 



where Ci is the molar concentration of an ion. The 
linear term, /3/, in eq (1) was introduced by Hiickel 
[12] to account for changes in the dielectric constant 
of the medium as a function of/. There is no accepted 
value for the parameter /3. In all of our previous solu- 
bility studies of calcium phosphates [1-7] the value 
of /3 was set equal to zero, which is appropriate for 
solutions having low ionic strengths (< 0.1). 

b. Solubility Product Constant 

The solubility product constant of CaHP0 4 • 2H 2 
is defined as 

K sp =(Ca*+)(HP02-) 

= [Ca 2 +] [HP0 2 4 -]y Ca 2+ • y H Po r (3) 

in which the quantities in parentheses represent molar 
ion activities and in brackets represent molar concen- 
trations. 
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Apparent solubility products, £ s ' p , for each solution 
were calculated initially from the raw data using a 
standard iterative procedure [13]; the needed activity 
coefficients were calculated from eq (1) with (3 = 0. 
Values for ionization constants for H3PO4 and H 2 0, 
the association constants for the ion pairs CaHPOJ 

and CaH 2 PO+, and the constants A and B in eq (1) 
were the same as those used by Gregory et al. [1]. 
The possible presence of the ion pair NaHPO" was 
neglected in the calculation of K' svt . Henceforth any 
other quantity calculated in this manner (/3 = and 
[NaHPO^^O) is indicated as a primed quantity. 
Final values of solubility products, K sp , for each 
solution were calculated similarly, except that statis- 
tically derived values for f3 in eq (1) and for Km , the 
stability constant of ion pair NaHPO" (see sec. 2.3. d), 

were incorporated into the computations. The stability 
constant for the ion pair NaHPOj is defined as 



Km- 



[NaHP0 4 ]--y NaHP04 _ 



(4) 



[Na + ][HPO r ]-y, Na+ -y Hpor 

Ion size parameters for Na+ and CI" were taken as 4.0 
and 3.0 [14]; the value used for NaHP0 4 -, 4.0, was 
the same as for H 2 P04~. Henceforth, quantities 
calculated taking into account the estimated values 
for j8 and Ks are denoted by symbols without primes. 

c. Electroneutrality Unbalance and Ionic Strength 

For the system Ca(OH) 2 — H3PO4— NaCl— H 2 0. 
electroneutrality unbalance is defined here, 



U(±)=-2CiZi 



(5) 



in which the sum includes the ions Ca 2 + ; H + ; CaH 2 
POj; Na+; NaHP0 4 ~; OH; CI"; H 2 P0 4 -; HPO*-; 
and PC>4 _ . Apparent and final values of the unbalances, 
U(±)' and U (±), as defined in eq (5), and of the ionic 
strengths, /' and /, as defined in eq (2), were also 
calculated. 

d. Statistical Evaluation of fi, K N , and K° p 

Values of the parameter ft in eq (1), the stability 
constant of ion pair NaHPO^, and an "overall" 
adjusted value of solubility product constant, Xj? p , 
for CaHP04 * 2H 2 were determined from the raw 
data utilizing a generalized least-squares approach 
based on the procedure given by Deming[15] and Guest 
[16]; the application of this method to solubility data 
of CaHPC>4 * 2H 2 in the ternary system has been 
described in detail by Gregory et al. [1]. 

In brief, the raw data were adjusted giving consid- 
eration to the weights of the measurements, so that 
the sum of squares of the weighted adjustments in the 
measurements was minimal. The weights were taken 
as inversely proportional to the squares of the standard 
errors in the measurements (Ca^ ± 2.5%; P and P* 
± 1.5%; NaCl ± 1.0%; pH ± 0.01). These adjust- 
ments were subjected to the condition that the ad- 



justed measurements and the final estimated param- 
eters would fulfill the following restrictions of the ionic 
model (in which the Debye-Huckel equation is used to 
calculate single-ion activity coefficients): 

1 ~2^ L < Zi (6) 

logy/" = {— AZ\ VF*/(1 + Bat V7*)} +/?/* (7) 



P*- 



. p* 



1.0 



(8) 



2[Ca 2+ ]* + [H + ]*+[CaH 2 P0 4 + ]*+[Na + ]* 

= [OH"l* + [NaHPOri* +[C1T + [H 2 P0 4 T 

+ 2[HP02-]* + 3[P03-]* (9) 

k:= [Ca*+j* -rHPon*-y*a 2+ • y™\- W 

[NaHPQ 4 -]*-yS.HPO < - 



Kn — 



[Na + ]*-[HP0 4 -]*- TH * ,_-y* a+ 



(ID 



In eqs (6-11), the starred quantities are either the 
final statistically adjusted measurements or the 
quantities derived from such values. Equations 
(8) and (9) are the conditions for dissolution stoichiom- 
etry for CaHP04*2H 2 and electroneutrality balance 
for the solutions, respectively. 

In a separate, but similar calculation, the Davies 
equation (eq 12) for obtaining ion-activity coefficients 
was used instead of the Debye-Huckel equation 
(eq7): 

log yf = { -AZ\ \f*/(l+ \7*)} + MZfI*. (12) 



In this equation A is an adjustable parameter. The 
previously accepted value, 0.2, for A [9] has been 
revised [17] to a new value, 0.3, based on newer 
compilations of activity data. Nevertheless, an inde- 
pendent estimate for its value was obtained from the 
present solubility data. 

3. Results 

The compositions of initial and equilibrated solu- 
tions, and the dissolution stoichiometry of the solid, 
Car/(Pf — P ), are given in table 1. The equilibrium pH 
values in sets A, B, C, and D increased with increasing 
amounts of NaCl; this is probably due to the effect 
of increasing ionic strength on the ion activity coeffi- 
cient of H + . This effect is masked in set E because 
formation of NaHPO^" has tendency to reduce pH; 
thus, the net effect of the two counter-balancing proc- 
esses is a near constancy in pH in solutions E-2 to 
E-6. Calcium and phosphorus concentrations were 
nearly constant within the sets A, B, and C; however, 
appreciable increases in the concentrations were 
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observed in the sets D and E with increasing amounts 
of NaCl. Stoichiometries of the dissolved solid were 
close to unity for all the sets (last column in table 1) 
except for small deviations in solutions C—l to C-5 
and D-6. 

Statistically derived values of parameters /3 and 
A in the Debye-Huckel and Davies equations, respec- 
tively, and the values of Kn and X s ° p obtained with the 
two models are given in table 2. 

Apparent and final values of ionic strength, electro- 
neutrality unbalance, and solubility product for 
CaHP0 4 * 2H 2 are given in table 3; these were calcu- 
lated with the use of the Debye-Huckel equation. 
Percentage of total phosphorus bound to Na + in the 



form of the ion pair NaHPOr is given in the last 
column of table 3. Only / values are listed, because 
/' never differed from / by more than 0.5 percent. 

Table 2. Statistically derived values of parameters in the ion 
activity coefficient equations, stability constant for 
NaHPOj and solubility product constant for CaHP0 4 - 
2H 2 



Derived 
quantity 


Debye-HiickeT Equation 


Davies Equation 


Equation 

parameter 
Kn 

K S U P x 10 7 


j8 = 0.0626 ±0.0054 

7.0±2.4 Lmol- 1 
2.48 ±0.02 mol 2 A~ 2 


A = 0.228 ±0.007 

9.4 ±5.2 Lmol- 1 
2.25 ±0.02 mol 2 .l" 2 



Table 3. Calculated values of electroneutrality unbalance, ionic strength, solubility 
product constant, and percent of total phosphorus bound to sodium ions 



Solution 

No. 



Ionic Model With Debye-Huckel Equation 



= 0.0; K N = 0.0 



U'(±) 
MxlO 3 



Kn X10 7 



P = 0.0626; K N =1.0 



1, M 



U(±) 
Mx 103 



K, D x lO 7 



100 x 
[NaHP0 4 -]/P* 







Pc 


= 14.8 Mx 10 3 






A-l 


0.361 


2.52 


0.0422 


0.361 


2.55 


0.0 


A-2 


.237 


2.42 


.0928 


.253 


2.48 


.051 


A-3 


.0936 


2.39 


.143 


.128 


2.49 


.11 


A-4 


-.0633 


2.36 


.194 


-.00765 


2.48 


.18 


A-5 


-.190 


2.31 


.294 


-.0870 


2.50 


.34 


A-6 


-.179 


2.11 


.545 


.0471 


2.44 


.74 



Po = 13.2Mxl03 



B-l 


-0.194 


2.50 


0.0378 


-0.193 


2.53 


0.0 


B-4 


-.514 


2.39 


.189 


-.452 


2.51 


.23 


B-6 


-.407 


2.23 


.541 


-.150 


2.57 


.93 







Po 


= 10.5Mxl0 3 






C-l 


0.670 


2.53 


0.0309 0.670 


2.56 


0.0 


C-2 


.653 


2.49 


.0815 .675 


2.55 


.094 


C-3 


.517 


2.48 


.132 


.565 


2.57 


.21 


C-4 


.263 


2.42 


.182 


.339 


2.53 


.34 


C-5 


.152 


2.39 


.283 


.291 


2.58 


.63 


C-6 


-.208 


2.20 


.534 


.0952 


2.52 


1.4 







Pp 


= 3.52Mxi0 3 






D-l 


0.0330 


2.44 


0.0110 


0.0332 


2.45 


0.0 


D-2 


.00260 


2.41 


.0613 


.0497 


2.44 


.65 


D-3 


-.0502 


2.38 


.112 


.0525 


2.42 


1.4 


D-4 


-.0912 


2.37 


.163 


.0705 


2.43 


2.2 


D-5 


-.224 


2.32 


.264 


.0540 


2.40 


3.8 


D-6 


-.0743 


2.30 


.515 


.511 


2.47 


7.6 







Po 


= 1.21MX10 3 






E-l 


-0.0488 


2.52 


0.00485 


-0.0488 


2.52 


0.0 


E-2 


-.192 


2.50 


.0564 


-.119 


2.47 


2.9 


E-3 


-.299 


2.57 


.107 


-.155 


2.51 


5.5 


E-4 


-.336 


2.55 


.158 


-.128 


2.46 


7.7 


E-5 


-.424 


2.56 


.259 


-.0976 


2.45 


11.6 


E-6 


-.676 


2.59 


.511 


-.0959 


2.44 


18.7 



Average A: SD = 2.49±0.05X 10-' mol 2 . 1" 2 
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4. Discussion 

Calculation of ion activity products for solutions 
of calcium phosphates with ionic strengths of about 
0.16 is of interest to investigators involved in the 
study of biological calcification. In dilute calcium 
phosphate systems (/<0.1), the Debye-Hiickel 
and the Davies expressions (/3 = 0.0 or A = 0.3) 
have given consistently reliable ion activity products 
as judged by the constancy of the solubility products 
in solutions saturated with respect to solid phases 
such as CaHP0 4 • 2H 2 [1], CaHP0 4 [2], /3-Ca 3 (P0 4 )2 
[3], and Ca 5 (P0 4 ) 3 OH [4-7]. The data presented here 
test the applicability of these expressions for solutions 
with ionic strengths in the range 0.005-0.55. These 
data also show that significant concentrations of the 
ion pair NaHPO^ form at high concentrations of 
Na + and at higher pH's; this obviously has considerable 
significance relative to the activities of phosphate 
ions in fluids such as blood. 

4.1. Solubility Product Constant 

The data presented in table 3 and figures la and 2a 
show that at higher values of / the Debye-Hiickel 




6.58 



~o — ar 



(3 = 0.062u 



(b) 



0.1 0.2 0.3 0.4 0.5 0.6 

IONIC STRENGTH, I 

Figure 1. (a) Plot of —log K' si) versus ionic strength; Debye-Huckel 
equation with ft = 0.0 was used and ion pair NaHPO^was omitted 
from the calculations. 

Key tor the symbols: © bet A; Q Set B; & Set C; 3 Set D; and • Set E ot table 3. 

(b) Plot of —log K sp versus ionic strength; the Debye-Huckel equa- 
tion with (3 = 0.0626 and the value 7.0 for the stability constant of 
the ion pair NAHPO4 were used in the calculations. 

The broken line corresponds to a mean —log K sp value of 6.603; the standard error in the 
—\ogK sp values (±0.009) is shown by the bar to the right of the points. 



equation (/3 = 0) and the Davies equation (A = 0.3) 
do not give constant pK' sp ; instead these products 
decrease (figure la) or first decrease and then increase 
(figure 2a) as a function of increasing ionic strength. 
We attribute these deviations (i) to failure of the Debye- 
Hiickel and Davies expressions with /3=0 and A=0.3, 
respectively, at high ionic strengths, and (ii) to forma- 
tion of the ion pair NaHP0 4 ~. 

The final adjustments required to fulfill the restric- 
tions of the ionic model (eqs 6-11) with Debye-Hiickel 
activity coefficients were, in general, within the 
standard errors of the measurements and were 
randomly positive and negative. However, for solu- 
tions C-l to C-5, the adjustments for phosphorus 
concentrations were all negative (—2.6, —2.6, 
— 2.4, — 1.6, and — 1.9 % respectively) and larger than 
the standard error, 1.5 percent. Solution D-6 required 
a large positive adjustment (3.8%) in the calcium 
measurement (standard error 2.5%). Slightly abnormal 
dissolution stoichiometrics in these particular solutions 
(table 1) are explainable on the basis of the possible 
systematic errors in the phosphorus concentrations 
for the solutions C-l to C-5 and in the calcium con- 
centration for the solution D-6. 

The results of the final adjustments obtained with 
the Davies expression (eqs (6) and (8-12)), although 
not listed, were similar to those described above and 
the same conclusions can be drawn from them. 

Plots of pK' sv (-log K' sv ) and pK sl) (-log K sp ) 
against ionic strength are given in figures 1 and 2 
for calculations made with the Debye-Hiickel and 
Davies equations, respectively. 

When the statistically estimated value (0.0626) for 
the parameter /3 is used in the linear term of the Debye- 
Hiickel equation and when the presence of NaHPO 4 
is incorporated into the calculations, the final values 
of pK sv show greatly improved constancy (fig. lb). 
Note especially that the trend in the points apparent 
in figure la, is eliminated; the scatter of the points of 
series E (filled circles) is reduced through the incor- 
poration of the ion pair NaHPO^ into the model. 

The positive value of the parameter /3 obtained in 
this study (table 2) is in accord with the view that an 
increase in / decreases the dielectric constant of the 
medium, as was shown for 14 other electrolyte systems 
by Hasted et al. [19]. Moreno et al. [10] reported a 
value (—0.082) for /3 in their solubility experiments 
with CaHP0 4 -2H 2 in H3PO4 solutions containing 
NaCl. However, only three solutions containing NaCl 
were used in their study; the maximum NaCl concen- 
tration was only 0.1045 M, and the presence of ion 
pairs between calcium and phosphate species was not 
taken into account. However, we recalculated their 
data, taking into account the presence of the ion pairs 
CaHPO!! and CaH 2 P0 4 + , with the finding that the 
K sp became essentially constant and /3 in the Debye- 
Hiickel expression became non-determinant. 

Robinson and Stokes [20] have used a value of 
/3 = 0.055 in the Debye-Hiickel equation to describe 
experimental values of activity coefficients of NaCl 
solutions up to 2 M; the value, 0.0626, obtained in 
this study compares reasonably well with theirs. 
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Our value (table 2) for KJ P , 2.48 ± 0.02 X 10" 7 , 
obtained with /3 = 0.0626 compares well with the 
value, 2.56 ± 0.05 X 10" 7 , obtained by Gregory et al 
[1] through use of a similar statistical treatment of 
data (except that /3 was set equal to zero) for the 
ternary system in which the ionic strengths were too 
small for the term fil to make a significant difference. 

Similar calculations of K' sp and K sv were also made 
with the Davies equation; these results are shown in 
figure 2. In this case, the pK' SJ> values decreased as 
a function of /, ■ up to / = 0.18 and then increased 
rapidly. The pK sf > values obtained with the Davies 
equation (fig. 2b), contrary to those in figure lb ob- 
tained with the Debye-Hiickel equation, show an ap- 
parent trend. For this reason, the Debye-Hiickel ex- 
pression with (5 = 0.0626 appears to be preferable to 
the Davies equation, at least for this system. It should 
be noted that the persisting trend in figure 2 may be 
due to the absence of the ion size parameters, at, in 
the Davies equation. 

4.2. Electroneutratity Unbalances 

Theoretically, the electroneutrality unbalances, 
U(±)' and U(±) should be zero. However, in prac- 
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FIGURE 2. (a) Plot of —log K' sp versus ionic strength; Davies equation 
with A =0.3 was used and the ion pair NaHPO-r were used in the 
calculations. 

Key for the symbols: Same as in figure 1. 

(6) Plot of —log K sp versus ionic strength; the Davies equation with 
A = 0.228 and a value of 9. 4 for the stability constant of the ion pair 
NaHPOj were used in the calculations. 

The broken line corresponds to a mean —log K sp value of 6.647; the standard error in the 
— log K aD values (±0.023) is shown by the bar to the right of the points. 



tice there are at least four reasons why the unbalances 
differ from zero. (1) Experimental errors in the pH, 
and concentrations of calcium and phosphorus are 
always present; (2) the presence of unknown ionic 
impurities are not taken into account in the calcula- 
tion of unbalances; (3) an ionic species such as 
NaHPOj may be omitted from the working model even 
though it may have a significant concentration; and 
(4) errors in the relative concentrations of the various 
phosphate ions may result from the use of incorrect 
dissociation constants or from deficiencies in the 
expressions for activity coefficients. 

The unbalances listed in table 3 are small and ran- 
dom in sign for the sets A and B. The apparent trend 
in U(±)' in the sets D and E — almost all values are 
negative and increasing with increasing amounts of 
ion pair formation — is attributed to the omission of 
the ion pair NaHPOj in calculation of the values of 
f/(±)'; the corresponding U(±) values, which take 
into account the ion pair NaHPO^, do not show this 
trend, and, thus support the existence of that ion pair. 
The abnormal value of the unbalance in solution D-6 
and a trend in the unbalance of solutions C~\ to C-5 
may be caused by experimental errors as mentioned 
in the section 3. The inclusion of the linear term 
(0.0626 /) in the Debye-Hiickel equation does not 
change the U(±) ' values significantly. 



4.3. Dissolution Stoichiometry 

The dissolution stoichiometry of CaHP0 4 -2H 2 0. 
Ca*/(P*— Po) should be 1.0. In most of the cases 
(table 1) this stoichiometry is close to 1.0 and within 
the standard error ± 0.03. 



4.4. Metasfability With Respect to Other Calcium 
Phosphates 

All the equilibrated solutions in this study were 
supersaturated with respect to Ca 5 (P0 4 ) 3 OH (solubility 
product constant = 4T8 X 10 -59 [5-7]); the activity 
products (Ca 2 +) 5 (POJ-) 8 (OH"), varied from 
2.7 X lO" 58 in A-l to 3.6 X 10 52 in E-6. The solu- 
tions A-l to D-6 were undersaturated, and the 
solutions E~~l to E-6 supersaturated, with respect to 
/3 — Ca 3 (P0 4 )2 (solubility product constant =1.4 
X lO" 29 [3]); the activity product, (Ca 2+ ) 3 (PO?-) 2 , 
varied from 5.5 X 10 32 to 6.1 X 10" 30 in A-l to D-6 
and 3.5 X 10" 29 to 6.31 X 10~ 29 in E~l to E-6. All the 
solutionswere undersaturated with Ca8H 2t (P04)6 , 5H 2 
(solubility product constant = 1.1 X 10" 94 [18]); the 
activity product, (Ca 2 +) 8 (H+) 2 (P0 3 -) 6 , varied from 
4.1 X 10- 101 to 4.9 X 10~ 95 in A-l to £--6. The above 
activity products were calculated with the Debye- 
Hiickel equation (p = 0.0626 and K N = 7.0) . 

Thus the solubility study reported here involvec 
metastable equilibria; this situation also existed ii 
previous solubility studies with CaHP0 4 *2H 2 [1] 
CaHP0 4 [2], and - Ca 3 (P0 4 ) 2 [3]. The effects ol 
crystallization of Ca 5 (P0 4 ) 3 OH or p - Ca 3 (P0 4 ) 2 
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during the solubility experiments are not important 
because of very slow growth rates of these salts. Sig- 
nificant precipitation of other calcium phosphate 
salt in the solution would cause deviations in the ap- 
parent dissolution stoichiometry of CaHPO^FbO. 
However, the values of U(±)' and U{±) would not be 
affected. Even though the greatest potential for precipi- 
tation of Ca 5 (PO^OH or /3-Ca 3 (P04)2 is in solutions 
E-l to E~6, close to ideal stoichiometry of dissolu- 
tion is observed in these solutions, confirming the 
earlier conclusion that the deviations in C~ 1 to 
C-5 and D— 6 must have been caused by factors other 
than precipitation of other calcium phosphate salts. 

4.5. Formation of the Ion Pair NaHPOj 

Formation of this ion pair should become significant 
at high Na + and HPOl" activities. Percentages of total 
phosphorus in the form of this ion pair are given in the 
last column of table 3. It is seen that a maximum of 
19 percent of the phosphorus is bound in the form of 
this ion pair (solution E-6). For comparison, a serum 
with Na+- 0.140 M, Ca, = 2.45 X 10~ 3 M, P,= 1.33 
X10 3 M, pH-^7.4 and 7=0.16 would have about 
18 percent of its phosphorus in the form NaHPOj • This, 
of course, would be a maximum value because of the 
presence of other ion pairs and complexes, but this cal- 
culation reveals that this ion pair should be given 
consideration. 

The stability constant obtained with the Davies equa- 
tion, 9.4 ±5.2, (table 2), for NaHPOj, is within the 
limits of error of that obtained with the Debye-Hiickel 
equation, 7.0 ±2.4. 

An "apparent" stability constant value, 4.0, was 
reported by Smith and Alberty [21] in a study carried 
out at constant ionic strength of 0.2. Through the use of 
this ionic strength and the Debye Hiickel equation 
(j8 = 0.0626), a corrected value 18 is obtained. This 
cannot be compared directly with our value 7.0 because 
its uncertainty is not known, but at least it has the same 
order of magnitude. 

Formation of this ion pair would increase the dissolu- 
tion of a calcium phosphate in a sodium-containing 
solution. This is apparent in the increases in the con- 
centrations of calcium and phosphorus with NaCl 
concentration in Sets D and E as compared to setsy4 to 
C (table 1). 

The greater formation of this ion pair in set £, as 
compared to sets A to C. is also manifested by the fact 
that the pH increases with increasing ionic strength 
much less slowly in set E than in the other sets. Here, 



formation of the ion pair NaHPOj has the effect of 
increasing the overall dissociation of the phosphoric 
acid. Thus, the presence of this ion pair may have a 
significant role in physiological processes. 



The authors thank E. M. S. Miletta for technical 
assistance. 
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